Rayleigh and Solberg criteria reversal near black holes: 
the optical geometry explanation 
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■^q , Summary 

The familiar Newtonian version of the Rayleigh criterion demands that for 
dynamical stability the specific angular momentum should increase with the 
increasing circumferential radius of circular trajectories of matter. However, 
sufficiently close to a black hole the Rayleigh criterion reverses: in stable flu- 
C~- ids, the angular momentum must decrease with the increasing circumferential 

radius. The geometrical reason for this reversal is that the space is so very 
strongly warped that it turns inside-out. 
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X : 1 INTRODUCTION 

Abramowicz and Prasanna (1990) found that, surprisingly, in the region* r < 3M around of a 
nonrotating (Schwarzschild) black hole with the mass M, the necessary condition for stability 
of axially symmetric, rotating isentropic and incompressible fluid against small, adiabatic and 
axially symmetric perturbations is that the specific (per unit mass) angular momentum of the 
fluid j should decrease with the increasing circumferential radius of the circular trajectories 
of the fluid r, 



dr 



<0. (1.1) 



The result was very surprising indeed, because it seemingly contradicted the old classical 
result of Rayleigh that, in Newton's theory, the necessary condition for stability should be, 

f)>°- (1-2) 
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Already a long time ago the Rayleigh criterion (1.2) was generalized to isentropic, compress- 
ible fluids by Solberg, and to non-isentropic fluids by H0iland. Rayleigh and Solberg criteria, 
which both have identical form (1.2), can be derived from the more general H0iland criterion 
as special cases (see Tassoul, 1978 for a detailed description of these criteria). Relativistic 
version of the H0iland criterion has been derived with a great care, and in a very mathemati- 
cally elegant way, in an important paper by Seguin (1975), but the effect of the Rayleigh and 
Solberg criteria reversal has not been noticed there. 

The reversal of Rayleigh and Solberg stability criteria provides yet another nontrivial example 
of an ambiguity of the geometrical meaning of "inside" and "outside" directions. I pointed 
out a few years ago (Abramowicz, 1992) that several counter intuitive effects in dynamics of 
bodies closely orbiting black holes could be understood as resulting from a strong warping 
of geometry that literally turns the space inside out! In the present paper, I show that the 
surprising effect of the Rayleigh and Solberg criteria reversal may also be explained the same 
way. 

2 THE LICHTENSTEIN THEOREM 

My discussion will be restricted to the motion on the equatorial symmetry plane of the 
rotating fluid. This restriction is not relevant for the results, but it greatly simplifies some 
arguments. 

Lichtenstein (1933) has proved that in Newton's theory stationary rotating fluid bodies always 
display a mirror symmetry with respect to reflections in a plane orthogonal to the rotating 
axis and containing the barycenter of the body (the equatorial plane) . For stationary, axially 
symmetric fluids, the equatorial plane is itself rotationally symmetric, and thus one may 
introduce on it a family of concentric circles around the symmetry centre. The circumferential 
radii r of these circles provide a convenient radial coordinate on the equatorial plane. 

In the early 1970s I was a Ph.D. student of Andrzej Trautman in Warsaw. He used to give his 
students a lot of freedom in performing research, and several of us have been indeed involved 
in doing all kind of esoteric theoretical physics in addition to our official Ph.D. subjects. 
Chandrasekhar who was often coming to visit Trautman in Warsaw, asked me once what 
I was working on in addition to my main subject. I answered, proudly, that I was trying 
to prove Lichtenstein' s theorem in general relativity. Chandra could not hide his obviously 
disapproving surprise, and only after a moment of silence replied very politely, "This sounds 
a bit too formal to me". Of course, when Chandra left, Trautman told me to stop with 
Lichtenstein and do something less formal. 

It is still unknown whether the Lichtenstein theorem is valid in Einstein's relativity, but 
all the known solutions of Einstein's field equations describing stationary, axially symmetric 
spacetimes do show this discrete symmetry. For example, the Kerr solution in Boyer-Lindquist 
coordinates is obviously symmetric with respect to reflection 9 — > it — 9 with respect to the 
equatorial plane, 9 = n/2, 



3 RAYLEIGH AND SOLBERG CRITERIA IN NEWTON'S THEORY 



In this Section I remind the familiar heuristic proof of the Rayleigh and Solberg criteria, first 
given by Randers (1942), and in the next Section I use a slightly modified proof to derive the 
relativistic version of the criteria. 

In a reference frame which corotates with the small fluid element under consideration, three 
forces are acting on the element: gravitational force G, centrifugal force Z, and applied force 
(e.g. pressure) T. On a particular circle r = r it is, 

G(ro) + Z(r ,Jo 2 )+r(ro) = 0. (3.1) 



Consider an adiabatic, axially symmetric perturbation in which the small element initially 
located at ro is displaced to r = r$ + dr. As the perturbation is adiabatic and axially 
symmetric, it does not change the specific angular momentum of the element. The (squared) 
angular momentum at the perturbed position r = r + dr is therefore Jq . On the other hand, 
the (squared) specific angular momentum needed at this very place for the balance of the 
forces (3.1) should be, 

J 2 = f(r + 5r) = J 2 + Sf = j 2 + (^j Sr. (3.2) 

Because j 2 ^ j 2 , there is unbalanced force ST at ro + Sr, 

T(r) + G(r) + Z(r 7 j 2 - 5j 2 ) = 5T. (3.3) 
Subtracting (3.1) from (3.3) and using (3.2) one gets, 

^-(flWV. (3.4) 



dp J \ dr 

Of course, the perturbation could be stable only if the unbalanced force brings the perturbed 
element back to its original position, i.e. only if 5T5r < 0. Thus, the necessary condition for 
stability yields, 

aZ ^^>0. (3.5) 



dj 2 J \ dr 

In the corotating frame, according to Newton's theory, Z(r, j 2 ) = j 2 /r 3 . Therefore, in 
Newton's theory, dZ/dj 2 is unconditionally positive, 

dZ\ 1 , s 

a ? j = ? >o. (3-6) 

and this is why the Rayleigh and Solberg criteria take the familiar form (1.2). However, 
Abramowicz and Prasanna (1990) have shown that, according to Einstein's theory of general 
relativity, it could be that dZ /dj 2 < 0. In this case the Rayleigh and Solberg criteria demand 
that for stability dj 2 /dr < 0. 



4 RAYLEIGH AND SOLBERG CRITERIA IN EINSTEIN'S THEORY 



The Rayleigh, Solberg and H0iland criteria are valid for stationary and axially symmetric 
equilibria of rotating fluids. Thus, their generalization to Einstein's theory should be derived 
in stationary and axially symmetric spacetimes. In such spacetimes two commuting Killing 
vector fields exist: asymptotically timelike Killing vector rf (corresponding to stationarity), 
and spacelike vector £ l with closed trajectories (corresponding to axial symmetry). The 
stationary observers are defined as those with velocities N l = 77* (1777) -1 ' 2 . More convenient 
are locally non rotating observers with four velocities, 



n , = e , V + aO, (4.1) 

e- 2 * = ( W )-u;(££), " = -M. (4.2) 

One may check that n l is a unit timelike vector which has zero vorticity and is hypersurfaces 
orthogonal. The hypersurfaces orthogonal to trajectories of n % are three dimensional spaces 
t = const with the metric, 

hik = 9ik - riiUk. (4.3) 

Because (n£) = 0, the three dimensional space is axially symmetric, and a circle in space can 
be defined as a trajectory of the projection of into space. A particle which moves along a 
circle in space has its four velocity in the spacetime given by, 

u i = 7 (n i + i;r i ). (4.4) 

Here r % = £ l /r, with r 2 = — (££) is a unit vector orthogonal to n\ and v is the orbital speed. 
From (uu) = (nn) = 1, (rr) = —1 one deduces that 7 = 1/(1 — -u 2 ) 1 / 2 , so that 7 is the 
Lorentz redshift factor. Note that the projection of the four velocity u 1 into the 3-D space 
gives, v = u l h k i = r yvr k and therefore v 2 = —(vv) = "f 2 v 2 . While the velocity v changes 
between ±1, the "velocity" v changes between ±00. 

Abramowicz, Nurowski and Wex (1993) demonstrated that in the corotating reference frame 
of matter that moves according to (4.4), the acceleration = u k VkUi equals, 

v 2 _ 
a, = V,$ + 7 2 — ViR-<y 2 vRV z uj, (4.5) 

where R is the radius of gyration R defined as (Abramowicz, Miller and Stuchlfk, 1993), 

# = re * = [-(££)] 1/2 e* (4.6) 

Abramowicz, Nurowski and Wex explained that the acceleration formula (4.5) suggests a very 
natural and convenient way to covariantly define gravitational, centrifugal and Lense-Thirring 
( "Coriollis" ) forces, 

v 2 

Qi = Vi^, Z l =~ f 2 —V l R, C i = - 1 2 vRV i oj. (4.7) 
ti 

In order to repeat dervivation of the stability criterion in a way similat to that discussed 
in the previous Section, one needs to know how to express the orbital speed by angular 



A 



momentum which is conserved during an axially symmetric, non-stationary perturbation. In 
relativity, exactly as in Newton's theory, orbital speed v, and specific angular momentum per 
unit energy £, 

'-8- (48) 

are connected by the "gyration equation" , 

v = L (4.9) 

The specific angular momentum per unit energy £ is not, in general, conserved by axially 
symmetric perturbations, because such perturbations need not to be stationary. The quantity 
conserved in axially symmetric perturbations (including nonstationary ones) is the specific 
angular momentum per unit mass, 

j = ££, (4.10) 

where £ is the specific per unit mass energy. In Newton's theory j = £ for any fluid, but in 
relativity j ^ £. However, in the special case of isentropic fluid considered here, the relativistic 
von Zeipel theorem (Abramowicz, 1974) guarantees that £ = £(£), j = j(£), and 

di 2 £ 2 

Here O is the orbital angular velocity measured by stationary observers. It is connected to 
the orbital angular velocity measured by locally non rotating observers O by, O = O — u, 
where 

n= v 4 = - = ^- (4.12) 

£ R R 2 y 1 

Using general expressions (4.7) for inertial forces, gyration equations (4.9), and assuming that 
the fluid is isentropic, I write on the equatorial plane, 

= ^, (4-13) 

Z(r£Hj))= £2{j) ±™ (4 14) 

1 ' [3)) R 2 -£ 2 (j)R* dr' 1 } 

C(r£(j)) = - i{j)R2 ^ (4 15) 

It should be obvious that the acceleration formula (4.9) may now be written on the equatorial 
plane in the form similar to the Newtonian equation (3.1), 

g(r ) + Z (r , £ 2 (jo)) + C (r , £(jo)) + T(r ) = 0. (4.16) 



The additional Lense-Thirring force C is not present in Newtonian case, and one must re- 
member that the quantity conserved in the perturbation is not £q but jo. Taking into account 
that dj 2 /d£ 2 > according to (4.11), one writes the general form of the criterion in the form, 



After a tedious calculation that utilises gyration equation one can rearrange this and write, 



dlnR l + v 2 du\ (dl 2 



dr 20 dr J \ dr 



> 0. (4.18) 



Seguin (1975) wrote the criterion in a different, but equivalent form, which follows from (4.17) 
also after tedious calculations. In the notation adopted here it yields, 

dl _ R2 dto\ fdl\ ^ 
dr dr ) \dr ) 

Although (4.18) and (4.19) are fully equivalent, it is convenient to discuss physical and geo- 
metrical meaning of the criterion using (4.18). 

5 REVERSAL OF THE RAYLEIGH AND SOLBERG CRITERIA 

The reversal of the Rayleigh and Solberg criteria occurs when the quantity 

S(r,,)^,f (5.1) 

changes its sign, i.e. at the circle Q(r,v) = 0. I shall explain the physical and geometrical 
meaning of Q = by considering precession of gyroscopes moving on circular orbits, and 
circular free motion of photons. 

The spin of an orbiting gyroscope does precess. The precession rate with respect to the 
vector t 1 tangent to the circle is given by the Fermi-Walker derivative, ui* = u^VjUi — 
{u 3 a,i — UiCi^Tj. On the equatorial plane, this vector has only one nonzero component which 
equals (Abramowicz, 1993), 



* 3 

u = 7 



1 2 )R^ UJ 
dr 2 dr 



(5.2) 



The factor 7 3 in the front of the bracket is due to special relativistic effects and represents 
the well-known Thomas precession. The first term in the bracket represents the geodesic (de 
Sitter) precession and the second term is the gravomagnetic (Lense-Thirring) precession 1 '" . In 
a weak gravitational field, dR/dr ~ 1 > 0, and dR/dr ^> \du>/dr\, so that u* and O have 
opposite signs. This is exactly what is expected by intuition which tells that an orbiting 
gyroscope must precess backward with respect to its orbital motion in order to always point 
to the same direction in space. Thus, gyroscopes orbiting clockwise precess anti-clockwise 
and vice versa. However, surprisingly, the bracket in (5.2) may change its sign. When this 
happens, the sense of the precession is exactly opposite to what is expected: gyroscopes 
are precessing forward — those on clockwise orbits are precessing clockwise, and those on 
anti-clockwise orbits are precessing anti-clockwise. (Rindler and Perlick, 1990; Abramowicz 



* There is a considerable confusion in the literature about how one should interpret precession of an orbiting 
gyroscope in terms of different kinds of simple precessions. In particular, some authors have previously been dis- 
cussing a rather unfortunate concept of the "gravitational Thomas precession" . This was because in approximate 
formulae used previously by several authors, the simple and unique division provided by (5.2) was not apparent. 



1990, 1993; Nayak and Vishveshwara 1998). It is not difficult to see that the reversal of the 
Rayleigh criterion and the reversal of the sense of gyroscope precession go hand in glove. 
Indeed, (5.2) may be written as, 

^* = -7 3 ^Q, (5-3) 

which means that both reversals are governed by exactly the same condition Q = 0. 

To see what does this condition really mean, one should consider the free (geodesic) motion 
of photons which is described by general equations a\ = 0, and v = ±1. It follows directly 
from (4.9), that they are equivalent to 

One should note that (5.3) is the same as the ultra-relativistic limit, v = ±1, of both (5.1) 
and (5.2). Thus, for the ultra-relativistic fluid, both reversals — of the Rayleigh and Solberg 
criteria and of the sense of the precession of gyroscopes - - occur at the location of the 
circular photon orbit. The conclusion that the reversal of the sense of gyroscope precession 
occurs at the location of a circular photon orbit has been previously reached by Nayak and 
Vishveshwara, 1998. 

6 DISCUSSION AND CONCLUSIONS 

Abramowicz, Carter and Lasota (1988) have noticed that dynamics of particles and photons 
looks simple when described in a conformally rescaled three dimensional geometry of space 
with the metric, 

hik = e~ 2<i> h ik . (6.1) 

The conformally rescaled metric (6.1) is called the optical reference geometry^. 
Consider a static spacetime, u> = 0. Its metric can be written in the form, 

ds 2 = e" 2 * (dt 2 + h ik dx l dx k ^ , (6.2) 

where the time coordinate t is invariantly defined as the synchronized time of the static 
observers, who have four velocities given by N l = n l . It is, = rf . Relativistic Fermat's 
principle assures that J dt on the light trajectory has an extremal value, 5 J dt = 0. This, 
together with dt 2 = —hikdx 1 dx k that follows from (6.2) and the fact that ds = for light, 
is equivalent to say that for any light trajectory x % = x l (p) in the space of optical reference 
geometry, 

i 

This means that light trajectories are geodesic lines in the optical reference geometry, and 
this is the reason for its name. 



" The optical reference geometry (6.1) has been introduced by Dowker and Kennedy (1978) and further discussed 
by Gibbons and Perry (1978), and Kennedy, Critchely, and Dowker (1980). It should be not confused with the optical 
geometry introduced by Trautman (1984) and Robinson and Trautman (1984) which is a totally different concept. 
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The circumference radius of a r = const circle equals, in the optical reference geometry, to the 
radius of gyration of the circle, because both radii are given by the same formula, R = re*. 
I shall now calculate the radius of the curvature of the circle in optical geometry, 1Z, using 
the well-known Frenet construction (see e.g. Synge and Schild, 1959), 

f l V,f fc = -A fc i. (6.4) 
K 

Here A is the unit, outside pointing, vector orthogonal to the circle, f l is the unit vector 
tangent to the circle, 

f* = e-*, fi = e*, (6.5) 
and the covariant derivative operator Vj obeys, 

r'V l f k = T l V,T k + V fc $ + r fc r l V,$. (6.6) 

From (6.4) I derive, 

OR R 



where e = sign(dR/dr). The formula for the gyroscope precession takes now the form, 

(6.8) 



* 3 

u> = 7 



In a curved space, consider a circle with circumferential radius R and curvature radius 71. If 
R^ Tl, a simple geometrical construction reveals that there is a nonzero deficit angle, 

^ = 2tt^^. (6.9) 

In the case of a very slow orbital motion, 7 = 1, in the static space, u(r) = 0, the precession 
of a gyroscope orbiting the circe with a steady angular velocity O = 2tt/T (T is the orbital 
period) is obviously given by, 

T 11 y J 

This is exactly the first term in formula (6.8). Thus, this term is indeed just exactly the 
geodesic precession — but "geodesic" in the sense of the optical reference geometry hik, not 
the directly projected geometry hik- 

In this article, I have described a rather counter intuitive effect of the reversal of the Rayleigh 
and Solberg stability criteria, and explained that it could be understood in terms of a very 
simple intuitive argument with a help of the optical reference geometry. In addition to this 
effect, the optical reference geometry has already successfully explained a score of similarly 
strange other effects. One should ask therefore the question: why is the optical reference 
geometry, and not the "real" geometry, so excidingley well fitted to behaviour of particles, 
fluids, charges, photons, gravitational waves, electromagnetic fields and gyroscopes? Are we 
missing something? Is the optical geometry just a particularly good map of the real curved 
geometry, or is the optical geometry itself the real geometry of space? 
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